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We study a simple open quantum system with a VT -symmetric defect potential as a prototype in 
order to illustrate a number of general features of PT -symmetric open quantum systems; however, 
the potential itself could be mimicked by a number of PT systems that have been experimentally 
studied quite recently. One key feature is the resonance in continuum (RIG), which appears in both 
the discrete spectrum and the scattering spectrum of such systems. The RIG wave function forms 
a standing wave extending throughout the spatial extent of the system, and in this sense represents 
a resonance between the open environment associated with the leads of our model and the central 
"PT-symmetric potential. We also illustrate that as one deforms the system parameters, the RIG 
may exit the continuum by splitting into a bound state and a virtual bound state at the band edge, a 
process which should be experimentally observable. We also study the exceptional points appearing 
in the discrete spectrum at which two eigenvalues coalesce; we categorize these as either EP2As, 
at which two real-valued solutions coalesce before becoming complex-valued, and EP2Bs, for which 
the two solutions are complex on either side of the exceptional point. The EP2As are associated 
with PT-symmetry breaking; we argue that these are more stable against parameter perturbation 
than the EP2Bs. We also study complex-valued solutions of the discrete spectrum for which the 
wave function is nevertheless spatially localized, something that is not allowed in traditional open 
quantum systems; we illustrate that these may form quasi-bound states in continuum (QBIGs) under 
some circumstances. We also study the scattering properties of the system, including states that 
support invisible propagation and some general features of perfect transmission states. We hnally 
use our model as a prototype for the construction of scattering states that satisfy PT-symmetric 
boundary conditions; while these states do not conserve the traditional probability current, we 
introduce the PT-current which is preserved. The perfect transmission states appear as a special 
case of the PT-symmetric scattering states. 


I. INTRODUCTION 

A. Two non-Hermitian systems: open quantum systems and PT-symmetric systems 

In the conventional formulation of quantum mechanics, the Hamiltonian operator H describing a given physical 
system is generally required to satisfy the Hermitian symmetry H = H\ a sufficient (but not necessary) condition to 
obtain a real-valued energy spectrum. Since the theory was originally developed, however, a number of researchers 
have found it useful to introduce non-Hermitian elements to the Hamiltonian, either as an extension of the original 
theory to accommodate certain physical situations HHia or as a useful reformulation in others In the latter 

case, various non-Hermitian Hamiltonians have been introduced to describe open quantum systems. 

Open quantum systems generally consist of a finite system coupled with an infinite environment, and thereby 
give rise to an energy spectrum with both discrete and continuous eigenvalues; the continuum is associated with the 
environmental degrees of freedom, while the discrete eigenvalues are a consequence of scattering due to the finite 
system. Some of the discrete eigenvalues can be complex, a signature of resonance phenomena in open systems. 
Resonances are associated with transient phenomena such as transport and exponential decay m HSl HDHIU [23 [571 
[3nj[3Il[Ml[35HS7] and may be viewed as generalized solutions of the Schrodinger equation with complex eigenvalues |331 
EiniEii] or as complex poles of the analytically continued S-matrix, among other perspectives |34l 138) . 

The reason why open quantum systems may accommodate complex eigenvalues can be summarized as follows. 
Eigenfunctions that are normalizable in open quantum systems, namely bound states and norm-preserving scatter¬ 
ing states, lie within the Hilbert space and can only have real eigenvalues. This corresponds to the fact that the 
Hamiltonian operator is Hermitian in the Hilbert space. However, even the standard Hamiltonian operator may be 
non-Hermitian in a space wider than the Hilbert space m- Open quantum systems indeed can harbor unnormaliz- 
able eigenfunctions, which lie outside the Hilbert space and can have complex eigenvalues depending on the boundary 
conditions. (Note, however, that we can still give a probabilistic interpretation for such eigenfunctions [371 155] 

While usually hidden in the boundary conditions, this non-Hermitian aspect of open quantum systems manifests 
itself when we trace out the continuous degrees of freedom associated with the environment; the resulting effective 
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Hamiltonian is then explicitly non-Hermitian. This effective Hamiltonian has only finite degrees of freedom remaining, 
corresponding to the discrete portion of the open quantum system, which is usually of primary interest. The first 
and most celebrated example in the literature may be the optical potential in nuclear physics. It was perhaps 
first introduced as a phenomenological potential but various researchers, Feshbach in particular, formulated it more 
rigorously [HlfT^HD] . 

In the case of the well-known tight-binding model this formulation leads to an energy-dependent effective potential; 
see Appendix C of Ref. [Ml. With the boundary condition of incoming energy we then have an effective Hamiltonian 
with a positive imaginary complex potential at the point where the discrete system couples to the environment, which 
represents an effective gain (or a source). On the other hand, a negative imaginary complex potential appears where 
the discrete system couples to the environment with the boundary condition of outgoing energy, which represents an 
effective loss (or a sink). 

As a recent development in the study of non-Hermitian physics, systems with both gain and loss have attracted 
a great deal of attention over the past two decades. Bender and Boettcher in 1998 demonstrated that one may 
relax Hermiticity in favor of T^T-symmetry (parity-time) in quantum mechanics and still obtain a real-valued energy 
spectrum in certain regions of parameter space mma. This has led some researchers to consider whether quantum 
mechanics could be reformulated in terms of VT symmetry; see, for example. Refs. I15H13 and particularly the 
references appearing in Refs. |481 149] . This theoretical question in turn inspired the idea of constructing physical 
systems that exhibit T^T-symmetry in the form of balanced gain and loss components arranged in a spatially-symmetric 
manner. 

A number of investigations have been carried out along these lines, both theoretically and experimentally, par¬ 
ticularly in the realm of optics |59llHT| . but also with examples in condensed matter physics jH], simple electronic 
circuits |63| . coupled mechanical oscillators |64H66] . and mesoscopic superconducting wires m- A number of intrigu¬ 
ing phenomena have been studied in the optical context, including power oscillations isniiniiia, double refraction HO] 
unidirectional invisibility |5E1 [5S11551 [5H] and localized states with novel transient behavior [23] ■ 

One central issue in the investigation of T^T-symmetric systems is VT symmetry-breaking. In many T^T-symmetric 
systems, one finds a transition between a phase in which all states are T^T-symmetric and a phase in which at least 
some states are not; the former is often referred to as the unbroken T^T-symmetric phase and the latter as the broken 
phase. At the T^T-symmetry breaking point, two real eigenvalues on the unbroken side coalesce and reappear on the 
other side of the transition as a complex-conjugate pair; their associated eigenfunctions are no longer T^T-symmetric 
individually, but only so as a pair (i.e. they appear as a state |'0) and its partner VT\ip)). We emphasize that at 
the transition point the eigenstates are not merely degenerate, but coalesce into a single state with a fixed universal 
phase between them Eaiin], as verified by experiment jZIHZS]. 

The T^T-symmetry breaking point, where the eigenstates coalesce, is an example of an exceptional point m- Similar 
transitions occur even in open quantum systems described by a Hamiltonian that is Hermitian within the Hilbert 
space. In this case the exceptional point is typically associated with the appearance of a resonance state (along with 
its anti-resonance partner) |7Klj77] after two real eigenvalues collide. While the large majority of studies on exceptional 
points appearing in the literature focus on the case of two coalescing eigenvalues (EP2s), the standard nomenclature 
is to refer to an exceptional point at which N eigenvalues coalescence as an EPA^ |I0J|7S|. In this paper we divide the 
EP2s into two further subcategories: we refer to an exceptional point at which two real-valued solutions meet to form 
complex conjugate partners as an EP2A; meanwhile we refer to an exceptional point at which two complex solutions 
with negative (positive) imaginary part coalesce to form two new solutions with negative (positive) imaginary part as 
an EP2B. 


B. "PT-symmetric open quantum system 

In this paper, we combine these two non-Hermitian systems in order to analyze a PT-symmetric open quantum 
system. Specifically, we incorporate a centralized VT -symmetric scattering potential iiT into an infinite tight-binding 
chain with otherwise real-valued site potentials. In the perspective given above, we can interpret this model as an 
otherwise standard open quantum system except that two sites are equipped with a direct environmental influence, 
one with -|-ir that injects energy into the chain, the other with —T that represents an energy drain. This may be 
realized as an optical lattice array in which one waveguide attenuates photon propagation (the ‘lossy’ component) and 
a second has a compensating amplifying character (the ‘gain’ component). We observe how the PT-symmetric gain 
and loss modify the usual open quantum system properties under two different boundary conditions: outgoing waves 
and scattering waves. For both of these we first consider the general case, including solutions that are PT-asymmetric, 
and then further investigate the solutions for which the boundary conditions themselves satisfy PT-symmetry. 

First we consider the boundary condition consisting of purely outgoing waves (often called the Siegert boundary 
condition) |27fE^ 15^ 1571 [75H55] . which yields the discrete spectrum for the system, including all bound states and 
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other solutions. We also observe the location of all exceptional points and other spectral features of interest. Here we 
demonstrate that for moderately small values of the T^T-parameter F, the spectral characteristics remain typical of 
traditional Hermitian open quantum systems. However, as we increase F explicitly non-Hermitian spectral properties 
emerge. 

We find a resonance state with vanishing decay width for certain specific values of F. In the context of a Hermitian 
open quantum system we would refer to this as a bound state in continuum (BIC) (see, e.g. Refs. |84f - IM) and 
references therein). While BICs typically appear owing to geometric effects and their wave functions discontinuously 
vanish outside a finite support, the present phenomenon results in a delocalized wave function with an eigenvalue 
that appears directly in the scattering continuum. For this reason, we refer to this state as a resonance in continuum 
(RIC). 

We further demonstrate the presence of localized states with complex eigenvalues that have recently been observed 
in an experiment and have since been considered in the theoretical works Refs. [smss]. We note that, unlike the 
RIC, these complex bound states appear over a wide range of parameter values, and, as observed in Ref. EH, the real 
part of the eigenvalue for these states may appear in the scattering continuum. Here we clarify that these localized 
states have complex conjugate values that sit in the first Riemann sheet in the complex energy plane, something that 
is not allowed in Hermitian open quantum systems. We also emphasize that while these states are indeed localized, 
they are not stationary states of the Hamiltonian. Instead, in an experiment they demonstrate either an amplifying 
or an absorbing characteristic |59| . However, given that the real part of these eigenvalues may reside within the 
continuum, in Ref. m the author classifies these states as a type of generalized BIC. By contrast, in this paper we 
emphasize that since these solutions are localized but non-stationary they would generally behave in a manner that 
is quite distinct from the usual concept of a BIC. That having been acknowledged, we further point out that there 
are some parameter ranges for which the imaginary part of the eigenvalues for these states will be very small, and 
hence they should take on a quasi-bound state behavior for these parameter values, similar to the quasi-bound state in 
continuum appearing in Ref. ESI mi. Specifically, these states should behave as bound states on time-scales t < F^/4, 
where the gain-loss defect parameter F exceeds the energy scale of the embedding optical bandwidth; we propose that 
these states might be detectable, for example, in a T^T-symmetric optical fiber loop array with a defect region [SS] 
that is modified to imitate our potential introduced in Sec. [ll| below (see Fig. [^. 

We then focus our attention on the ordinary bound state solutions appearing in our system and demonstrate that 
the wave function for these states satisfies T^T-symmetric boundary conditions. Further, we clarify that the wave 
function for virtual bound states (with real eigenvalue) is also T^T-symmetric, despite the fact that these states do 
not appear in the usual diagonalization scheme. 

We then consider the case of scattering wave boundary conditions. In the general case (T^T-asymmetric scattering 
waves) we observe that the parameter choices associated with the RIC result in a divergence in the reflection and 
transmission coefficients. This phenomenon has previously appeared in the literature in which it is referred to as a 
spectral singularity |92flM] and physically can be associated with both lasing and coherent perfect absorption |99H101| . 
We then demonstrate that a subset of the scattering wave solutions yield perfect transmission through the scattering 
region. In the special case in which the scattering potential is pure imaginary, we show that one can obtain perfect 
transmission for any continuum scattering states by appropriately choosing the value of F; this property approximately 
holds when small real-valued defects are introduced. We further demonstrate in this case that invisibility (perfect 
transmission with no scattering phase shift) can be obtained at discrete values within the continuum. 

In Sec.|T^below we present our prototype model for an open quantum system with a 'PT-symmetric defect potential. 
Then in Sec. |HI| we study the model under the boundary condition of outgoing waves, which yields the discrete 
spectrum associated with the defect potential. For the simplest case of a purely complex defect potential, we locate 
all exceptional points in the spectrum and characterize the properties of the spectrum in their vicinity; we further 
locate the RIC eigenvalues and write the associated wave function as an outgoing plane wave from the defect region. 
We also identify the parameter ranges that give rise to the localized states with complex eigenvalues and point out the 
situation in which some of these solutions might behave as quasi-bound states. In Sec. [Ill D] we generalize this picture 
by considering a potential with both real and imaginary defects. Here we demonstrate that as one deforms the system 
parameters, the RIC may exit the continuum by splitting into a bound state and a virtual bound state at the band 
edge; we believe that this point has not previously appeared in the literature. We note that traditional real-valued 
bound states also may appear for this more general potential. We study in closer detail the formal properties of the 
bound states in Sec. |IV[ demonstrating that they satisfy T^T-symmetric boundary conditions as expected. We also 
consider the CVT norm for these states, which we believe has only previously been investigated in closed VT systems. 

We then turn to the scattering boundary conditions in Sec.|V] which we use to characterize the RIC in greater detail. 
We also show that a subset of the scattering wave solutions give rise to perfect transmission through the scattering 
region, and in the case of a purely imaginary defect potential, there are two scattering solutions that support invisible 
signal propagation. We further demonstrate a connection between the localization transition in the discrete spectrum 
and the perfect transmission states that might be useful from the perspective of designing systems with predictable 
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Figure 1: Geometry for "PT-symmetric optical lattice with scattering potential given in Eq. ||^ 


transport properties. We also point out a possible application in the form of a ‘switch’ that is sensitive to invisible 
transmission originating from the left (right), but ignores such transmission from the right (left). In Sec. 


VI 


we 


demonstrate that a scattering wave solution can be obtained that itself satisfies T^T-symmetric boundary conditions. 
We also introduce the T^T-current, which is conserved for the (general) scattering wave solutions in our system, and 
which experiences a divergence associated with the perfect transmission states. We summarize our work and make 
concluding remarks in Sec.|VII| We also present some details of the calculations from the main text in two appendices. 


II. PT-SYMMETRIC OPTICAL LATTICE MODEL 

In the present paper, we study a tight-binding model with a VT -symmetric scattering defect potential, which can be 
realized as an optical lattice array or could be approximated by a modified version of the VT -symmetric optical fiber 
loop array with a defect studied in Ref. jH] or other systems appearing in the literature |bOI ll[)2| . Our tight-binding 
model takes the form 

OO 

H = - (|a: +l)(a;|-k |x)(x-k 1|)V(a;)|a;)(a;|, (1) 

X — — OC1 X 

in which the defect potential is specified as 

{ El -k jF for X = — 1, 

Eo for X = 0, (2) 

El — ir for X = — 1, 

where Eq, Ei, and F are all real, with V{x) = 0 otherwise, such that our scattering potential is confined to the central 
sites |x| < 1. The positive imaginary part of the complex potential contributes a factor exp[—i(f|F|)t] = exp(|F|t) to 
the time evolution, and hence is interpreted as being influenced by a particle bath that constantly injects energy as 
-kf|F|. That with a negative imaginary part is similarly interpreted as a particle bath that constantly drains energy 
as —*|F|. 

The off-diagonal part of the Hamiltonian Q is Hermitian (real symmetric) while the diagonal potential is not. 
It nonetheless satisfies the condition V{x)* = V{—x), which guarantees the system is TT-symmetric |51l 110311104] . 
Stated explicitly, the parity transformation V swaps the potentials at x = — 1 and x = -kl while the time reversal 
operator T (which is complex conjugation) flips them back to the original configuration. We note that several studies 
on TT-symmetric tight-binding models may be found in the literature, some of which are related to our model 

above EnunsHinH]. 

The Schrodinger equation = E\il)) for the Hamiltonian 0 can be written explicitly in the following way. 

First, let us consider the projection (x|7F|'0) = if(x|'(/') for the system component outside of the scattering potential 
|x| > 2, in which case V{x) = 0. We thus obtain 

—il}{x — 1) — i^ix -k 1) = ifi/)(x), (3) 

where '!/’(*) = Tho solution is given by ^^{x) = with the eigenvalue 

E{k) = —2cosk, (4) 

which defines the scattering continuum for our system in the range \E{k)\ < 2 with k G [—tt, tt]. To solve the eigenvalue 
problem in the scattering region, we hold the continuum dispersion E{k) and evaluate the Schrodinger equation for 
X = 0 and ±1, by which we obtain 


tp{-2) - ip{0) -k (ei -k *F)-!/'(-1) = E{k)jp{-1), 
-ipi-l) - -ipil) + Eo^/>(0) = E{k)'>p{0), 
-■0(2) - 0(0) -k (ei - fF)0(l) = E{k)tp{l). 


(5) 

( 6 ) 
(7) 











5 


first 




-2 


bound 


scatt 5 : 
conti ■ 


Le 


+2 


ring bound 
luum 


second 


anti-re 

X 

-2 


sonant 

X 


virtual branc 
(anti-bound) ^ 


reso 


+2 


Lb 


-K* 


^ virtual 
(anti-bound) 


nant 


(a) 


O 


SC 


atterin^ 


k 


continuum 


X 
X 

inti-resonant 


X 

a resonant 

,g : 


+7t 


B: 


(b) 


Figure 2: A typical distribution of the eigenvalues of the Hermitian tight-binding model on (a) a complex energy plane and (b) 
a complex k plane. 


A given solution ^/’(x) must satisfy these equations, subject to a specific choice for the boundary conditions. In Sec. Ill 
below we consider the boundary condition for resonant states that consist of purely outgoing waves, while in Sec. |V 
we consider the boundary conditions for scattering states. 

For later reference, let us present in Fig. [^a typical distribution of the eigenvalues of the Hermitian tight-binding 
model, that is, for F = 0. We will heavily make use of the terms in the figure. 

The complex E plane in Fig. |^a) consists of two Riemann sheets; they are connected by a branch cut that extends 
over the range —2 < E < +2. The first sheet corresponds to the upper half of the complex k plane in Fig. [^b), 
while the second sheet to the lower half. More specifically, the upper half (Imif > 0) of the first Riemann sheet 
in the complex E plane corresponds to the first quadrant of the complex k plane, the lower half (Im E < 0) to the 
second quadrant, the upper half of the second sheet to the third quadrant, and the lower half to the fourth quadrant. 
Notice that we can go, for example, from the upper half of the first sheet over to the lower half of the second sheet 
continuously through the branch cut, which corresponds to moving from the first quadrant to the fourth quadrant in 
the complex k plane. 

The scattering states continuously surround the branch cut on the real axis of the complex E plane. We hereafter 
refer to the scattering continuum as the energy band and to the end points of the continuum as the band edges, 
following the custom of condensed-matter physics. In the complex k plane, the scattering continuum is on the real 
axis, which is restricted to the first Brillouin zone — tt < fc < -I-tt; note that the line Refc = — tt is identified with the 
line Re fc = +tt as a result of the lattice periodicity. 

Bound states can exist on the first Riemann sheet below and above the energy band, that is, to the left {E < —2) 
and to the right {E > 2) of the scattering continuum. Those below the band lie on the positive imaginary axis of 
the complex fc plane, while those above the band lie on the positive part of the line Re fc = -I-tt. Notice that the 
bound states here have purely real eigenvalues; in other words, they never exist on the first and second quadrants of 
the complex fc plane except on the lines Re fc = 0 and Re fc = -I-tt. We will show below that once we introduce the 
non-Hermiticity, complex eigenvalues can appear in the first Riemann sheet (and therefore on the upper half of the 
complex fc plane); this is one critical difference between Hermitian and non-Hermitian open systems. 

The resonant states appear in the lower half of the second Riemann sheet, which is the fourth quadrant of the 
complex fc plane, while their anti-resonant partners reside on the upper half of the second sheet, which is the third 
quadrant of the complex fc plane. These are related to one another through time-reversal symmetry m- Virtual 
(or anti-bound) states can also appear to the left and right of the branch cut on the second Riemann sheet, which 
respectively correspond to the negative (Im fc < 0) parts of the lines Re fc = 0 and Re fc = -I-tt. 
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III. OUTGOING WAVES BOUNDARY GONDITION AND DISCRETE SPECTRUM 

In the present section our first consideration is the resonant states; there are several ways of computing these 
states |32]. We here use the Siegert boundary condition |27H291 EH |37l I79H83], which dictates that the system has 
outgoing waves only; this is equivalent to looking for all poles of the S matrix. The solutions of the resulting polynomial 
equation give the discrete eigenvalues associated with the scattering region, as shown below. 

Our purely outgoing wave function takes the form 


obtain 


-1 

0 


1 

r for a- < 

tpix) = < 

V'(o) 

for a; = 0, 

1 


for X > 1. 

r—1 

which 

brings Eqs. ^ 

-1 

0 

\ /V'{-1)\ 

^0 

-1 

V'(o) = 

-1 -A + 

£i — iT 

/ V V'(i) / 


( 8 ) 


^(0) 

^( 1 ) 


(9) 


in which we have introduced A = for convenience. In this notation the continuum dispersion Eq. takes the 
form 


E(A) =-(A + A-i). (10) 

We obtain non-trivial solutions for the discrete eigenvalues \j when the determinant of the matrix in Eq. ([^ 
vanishes. This is equivalent to the solutions of the quartic equation P{Xj) = 0 with 

P(A) = {el + r2) A^ + eo (e? + T") A^ - (l - e? - 2eoei - T^) A^ + (eg + 2ei) A + 1. (11) 


For a given solution Xj, the physical energy eigenvalue is determined from E{Xj) and the associated wave number 
is given as kj = —ilogXj. We emphasize that the number of solutions (four), is greater than the matrix dimension 
(three); this is because the matrix itself depends on the energy eigenvalue through the variable A. 

In the remainder of this Section we investigate the four discrete eigenvalue solutions of the quartic equation P{Xj) = 
0 in detail; first, we study the case of a purely imaginary defect potential with £o = £i = 0 in Sec. |III A| Here we 
locate all EPs and characterize the behavior of the spectrum in the vicinity of these points. We further identify the 
RIG and write the wave function of this state as a plane wave originating from the impurity sites. We also discuss 
the complex-valued localized states and their asymptotic localization properties as well as drawing attention to the 
parameter ranges in which some of these states will behave as quasi-bound states. 

In Sec. HID we generalize this picture to consider the case ei ^ 0 to illustrate two points (we keep eg = 0 for 
now). First we note that the EP2As still appear in the spectrum for the ei yf 0 case, while the EP2Bs vanish. 
This suggests that the EP2As may be more robust against parameter perturbations than the EP2Bs, on which we 
comment in relation to experimental results. Second, we demonstrate that as we increase the value of ei, one of the 
RICs approaches the band edge and eventually exits the continuum by splitting into a bound state and a virtual 
bound state. 


A. Discrete spectrum for eg = £i = 0: exceptional points (EPs), resonant states in continuum (RIGs) and 

quasi-bound states in continuum (QBICs) 


We first consider the discrete eigenvalue spectrum for the simplest case of our Hamiltonian Sq = Si = 0, for which 
the only non-homogeneous element remaining in the s yste m is the gain/loss pair governed by the VP parameter T. 
In this case, the quartic polynomial P(r) given in Eq. (11) simplifies to a quadratic in A^, yielding the four solutions 


Ai,4 

^ 2,3 


^V2r 


\/i - r2 -h Vi - 6r2 -t r4, 

- r2 - a/i - 6r2 -t 


( 12 ) 
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Figure 3: Discrete eigenvalue spectrum for the simplest case £o = £i = 0: (a) KeEj and (b) Im Ej against the PT-parameter 
F; (c) Refcj and (b) Im kj against F; parametric plots of (e) (Re (F), Im iJj (F)) and (f) (Re fej(F), Im A:j(r)) in the complex 
plane. In (e) and (f), the solid circles indicate some of the eigenvalues at F = 0, while the open circles indicate those in the 
limit F —>■ (xj; the arrows indicate how the eigenvalues evolve as F is increased from 0 to oo. 






















































We plot the real and imaginary parts of the resulting energy eigenvalues Ej = —Xj — A^- ^ as a function of F in Fig. 3 
(a) and (b), as well as the real and imaginary parts of the associated wave number kj = —ilogA^ in Fig. (c) ana 
(d). Figure ^(e) and (f) are parametric plots in the complex E plane and the complex k plane in the range of F > 0. 
In Fig. (^, (b) and (e), solutions plotted with full curves appear in the first Riemann sheet of the complex E 
plane while those plotted with a dotted curve appear in the second sheet; the former are the solutions with positive 
imaginary parts of kj and the latter are those with negative parts in Fig.[^(d) and (f). The first and second Riemann 
sheets of the complex E plane respectively corresponds to the upper- and lower halves of the complex k plane; a 
branch cut running from E = —2 to E = 2 connects the two Riemann sheets. We realize from the Siegert boundary 
condition (|^ that every solution on the first Riemann sheet has a positive imaginary part of the wave number and 
hence its wave function is bounded in x space, while every solution on the second Riemann sheet has a wave function 
that diverges along the leads of the optical array. 

We immediately observe one critical difference between Hermitian and T^T-symmetric open quantum systems: in 
the T^T-symmetric case, solutions with complex eigenvalues are allowed to appear in the first Riemann sheet, with 
localized wave functions. This is in stark contrast to the Hermitian case, in which complex-valued solutions are allowed 
to appear only in the second sheet, where they give rise to delocalized resonance and anti-resonance states. 

Let us summarize the evolution of the discrete eigenvalues from F = 0 to -|-oo along the lines of Fig. (e); the 
change for negative F is symmetric as this just amounts to swapping the gain and loss elements. At F = 0, one 
eigenvalue is at the lower edge of the continuum E = —2 and another at the upper edge E = +2 (solid circles in 
Fig. i (e) and (f)). There are also two eigenvalues at E = —oo and at E = -|-oo both on the real axis of the second 
Riemann sheet. 

As we increase F from 0, the eigenvalues at E = ±2 separate off from the band edges and move outward, while 
the eigenvalues at E = ±oo move inwards, all four along the real axis of the second Riemann sheet of the complex 
energy plane. These eigenstates are referred to as virtual bound states or anti-bound states in the sense that they are 
real-valued solutions that are spatially delocalized 133 EZl [ini]. The positive pair of solutions and the negative pair 
each coalesce at a point on the real axis of the second Riemann sheet at F = Fa = -s/S — 1, which is a second-order 
exceptional point. We label this point F = Fa as an EP2A and the region up until this point 0 < F < Fa as Region 
I. After passing the EP2A, all four eigenvalues become complex on the second sheet, forming two resonance/anti- 
resonance pairs symmetrically on the positive and negative sides. In the vicinity of the EP2A, the eigenvalues can be 
expanded in the characteristic form UlITT] 


El{r) = - 


^2(l + V2^ 


± 


2^/^V-l + V2 

for the resonance/anti-resonance pair with negative real part and 


Fi-F2 


E+(F) = W2 1 + V2 ± 


2i/4A/-l-Hy2 


^fi-F^ 


(13) 


(14) 


for the resonance/anti-resonance pair with positive real part. The derivation of these expressions is detailed in App.[A| 

We can regard Region I as the PT-unbroken phase and the EP2A at F = Fa as the PT-symmetry breaking point. 
As we see in Fig. Region I is the only continuous parameter region in which all discrete energy eigenvalues are real. 

As we continue to increase F, the complex eigenvalues eventually turn around and then return to the real energy 
axis at F = F^^j^ = 1. Although each pair of the energy eigenvalues are degenerate when they reach the real 
axis, their wave numbers are all distinct as can be seen in Fig. (c) and (f), and therefore this point represents a 
degeneracy in the standard sense, not a coalescence in the sense of the exceptional point. We refer to these states 
as resonances in continuum (RICs) for reasons described below in Sec. IIIB (also see Sec. VBl, and we refer to the 
region f a < F < FJ^jq as Region II. 

As we further increase F such that F > the four solutions pass through the branch cut running from E = —2 

to if = 2 and emerge on the first Riemann sheet of the E plane. This is equivalent to the observation that these 
solutions now have an effective wave number kj with positive imaginary part as shown in Fig. |^d) and (f). This 
implies that the wave function for these states 'tpji.E) ~ is localized, although the real part of the eigenvalues lie 

within the range —2 < Reif < 2 (see Fig. [^a)). These type of states were recently observed in an experiment based 
on light transmission through an effective T^T-symmetric array of optical fiber loops |59| in which they gave rise to a 
pair of exponentially growing and decaying localized states within the continuum. 

As F reaches the value F = Fb = 1 -I- V^, these states coalesce on the imaginary axis of the complex energy plane 
atF = fB = l-l- \/2, two on the positive side and the other two on the negative side, which is another second-order 
exceptional point (this time occurring in the first Riemann sheet). We refer to this point as an EP2B, because it 
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involves a pair of complex eigenvalues coalescing before becoming another pair of complex eigenvalues; we also refer 
to the region < F < Fb as Region IIL In the vicinity of the EP2B, the eigenvalues can be expanded as 


ifl 


-1 +V2 ± 


En{r) = - 


for the two eigenvalues with negative imaginary part, and 


,_./F2 - F2 

21/4vTT72'^ 


E+(T)=z^2 [-I + V 2 ) 


± 


21 / 4^1 + 72 






(15) 


(16) 


for the two with positive imaginary part, similar to the expressions near the EP2A above (see App. 0. 

After surpassing the EP2B, two eigenvalues move to the origin while the other two go off to ±zc», all on the 
imaginary axis of the first Riemann sheet of the complex energy plane. We refer to this region F > Fb as Region 
IV. Since F 3> 1 generally holds here, we may expand the solutions in Eqs. (12) in powers of 1/F to show that two 
of these solutions behave as ifi _4 « ±i(F — 2/F); note that in the limit F —)■ 00 , these two solutions asymptotically 
approach the simple value of the gain or loss component of the VT parameter F, as is indicated in Fig. [^b), where 
these two solutions (blue curves) approach the two diagonal (red) lines. Indeed, as shown in Appendix]^ the solution 
El ^ +iF is localized at site x = —1, while the solution E 4 ^ —zF is localized at the x = 1; hence, these two solutions 
gradually begin to mimic the original uncoupled gain/loss pair for large F. We comment further on the asymptotic 
localization properties of these states and show that the solutions E 2,3 behave as quasi-bound states in the continuum 
in Sec. IIIICI 

We emphasize that the physics in Regions I and II could arise in Hermitian open quantum systems as well; explicitly 
non-Hermitian properties appear in Regions III and IV with the appearance of the RIC and then the complex 
eigenvalues on the first Riemann sheet. 


B. Resonant state in continuum (RIC) 


Here we describe the resonant states in continuum (RICs) at the point F = F^jq in greater detail. As summarized 
above, the eigenvalues here appear on the real axis, embedded in the energy continuum that spans —2 < E < 2. At 
a glance, these states appear similar to bound states in continuum (BICs), which in Hermitian systems appear as 
resonances with vanishing decay width IMHM] . However, closer inspection reveals that these states are fundamentally 
different from BICs. 

For example, in the (Hermitian) double impurity open quantum system model studied in Ref. |88| it is shown that 
BICs appear as localized states between the two impurities; due to interference, the wave function for the BIC states 
vanishes identically outside of the impurity region. More generally, BICs often appear for geometrical reasons and 
hence are strictly confined in some spatial area. 

The present RICs, however, take the form 


i’mcix) = S for a; = 0 

1 gT*^(^-l)/4 for a; > 1 


(17) 


at F = F^jq, for the respective eigenvalues Eric = (with fcRic = ±37r/4) and Eric = — (with fcRic = ±7r/4). 

We refer to these points as resonant states in continuum (RICs) in part because the wave function for these states is 
delocalized as demonstrated in Eq. and because these states satisfy the Siegert boundary condition for outgoing 
waves. We will comment further on this naming convention in Sec. IVBl from the perspective of the scattering wave 
boundary conditions. We note that these are also equivalent to the spectral singularities that have previously appeared 
in the literature |92fI5H] . In Sec. HID we will also show that for the case ei yf 0, an RIC may approach the continuum 
edge and split into a bound state and a virtual bound state. However, we add one further brief comment here to 
emphasize that the RIC is not an exceptional point, as the eigenstates do not coalesce, having different wave numbers, 
and hence no fractional power expansion such as Eqs. (13 14) is possible in this case. 
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C. Quasi-bound states in continnum (QBICs) 

The solutions i?i ^4 from Region IV (or either pair of solutions from Region III) correspond to the localized states 
with complex eigenvalues that were recently experimentally observed in Ref. [59], in which the authors investigated 
light transmission through an effective T^T-symmetric optical lattice realized by periodically switching gain and loss 
in two optical fiber loops imiss]. As reported in Ref. j^, when a localized defect is introduced into the effective 
array (both a shift in VT pairing strength as well as a phase defect), a pair of localized complex conjugate modes 
appear within the continuum exhibiting exponential growth and decay in the power spectrum. Indeed, our solutions 
if 1,4 « ±*(r — 2/r) in Region IV appear directly in the center of the energy continuum (with Re if 1,4 = 0) and would 
also give rise to an exponential power output (growth or loss) as ~ 

While the author of Ref. [61] interprets these type of localized states with complex eigenvalue as examples of an 
effective BIC based on the fact that the real part of each solution may reside within the continuum, we note that since 
these states decay or grow exponentially, they would generally behave in a manner that is quite distinct from the 
usual concept of the BIC. However, the other pair of solutions if 2,3 also have the real part of the eigenvalue residing 
within the continuum, yet behave quite differently in Region IV. Indeed we can show that the eigenvalues for these two 
solutions behave as if 2,3 « ±i 2 /r^ such that the complex part of the eigenvalue for these states becomes arbitrarily 
small for increasing values of T. Hence, these states should behave as effective BICs on time scales satisfying t < r^/4, 
similar in concept to the quasi-bound state in continuum (QBIC) introduced in Refs. |9()L I91| . which are resonance 
states in the continuum with extremely long lifetime (also see Refs. [IIOIllll] '). 

As shown in App. the respective wave functions for the solutions £ 2,3 are exponentially localized around the 
site X = 0, while those for the solutions ifi ,4 are localized around the PT impurities at a; = ±1; we also show that 
the localization for the solutions if 1,4 is very narrow as it scales for T ^ 1 as l/logT, while that for the quasi-bound 
states if 2,3 is very broad, scaling as r^/2. We believe that these quasi-bound states should be observable, for example, 
in an experiment similar either to Ref. |59| or Ref. m in which the PT-symmetric defect potential is modified to 
mimic our potential appearing in Fig. 


D. Discrete spectrum for ei V 0: EP stability and RIG splitting at localization threshold 


As we relax the restriction Si = 0, most of the basic features that we observed in the simplest case in Sec. |HI A| -C 
remain, although these become somewhat distorted as shown for ei = 0.2 in Fig. |^a)-(d). Here we observe that the 
EP2As split into two pairs, one pair of which moves outwards and away from the origin on the F-axis while the other 
pair moves inwards towards the origin (for larger values of ei, the latter pair will eventually collide at the origin, 
before becoming complex-valued). While we obtained compact analytic expressions for the eigenvalue expansions in the 
vicinity of the EP2As for the case Si = 0, here those expressions become significantly more cumbersome. Nevertheless, 
following an intuitive generalization of the methods presented App. one may still easily obtain numerical versions 
of Eqs. (141 and in the vicinity of the EP2As in the more general case. 

On the other hand, the EP2Bs that we studied in Sec. |HI A| immediately vanish from the spectrum for ei ^ 0, as 
can be seen in Fig. gb) and the inset of Fig. ga); we can also see these coalescences vanish by comparing Fig. ge) 
and (f) (for ei = 0) with Fig. gc) and (d). Indeed, it can be shown that the EP2As also survive the generalization 
for £q 7 ^ 0, while the EP2Bs do not re-emerge. This seems to indicate that type EP2A exceptional points are more 
stable against parameter perturbations than those of type EP2B. We note that several experimental studies have 
been conducted in which an EP2A has been observed by simply passing directly through the exceptional point while 
varying a single parameter |53[ 1501 [03], but experimental observation of EP2Bs have tended to rely on encircling 
the exceptional point imiii] or mapping out the complex eigenvalue structure around the exceptional point in a 
two-dimensional parameter space |112j (although Ref. [72] provides an exception where the EP2B is observed more 
directly). Theoretically, we believe that the underlying reason for this is that the EP2As seem to vanish from the 
real parameter space only when they collide with another EP2A (see Ref. for another simple example where this 
occurs), while the EP2Bs do not need to collide with another EP in order to exit into the complex parameter space. 

The RICs meanwhile are also split apart in the parameter space, appearing at iFj^jp and iFj^jQ, given by 




i±kil 



(18) 


which we explicitly indicate by red crosses in Fig.ga) and (b). As we increase £i from 0, the RIG wave numbers for 
F^iq, which we refer to as in Fig. gd), move away from the ei = 0 values of 7 r /4 and 37 r /4 (previously shown in 
Fig. gf)) and approach 7 r /2 and the upper band edge tt, respectively. In the latter case, the RIG at Fj^jp eventually 
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Figure 4: Discrete eigenvalue spectrum for the case eo = 0. (a) Im Ej and (b) Imkj against the 'PT-parameter F for ei = 0.2; 
parametric plots of (c) (Rei?j(r), Imi5j(F)) and (d) (Re fcj (F), Imfcj(F)) in the complex plane for ei = 0.2; (e) ImEj and (f) 
Imkj vs. the "PT-parameter F for ei = 0.6. In (c) and (d), the solid circles indicate some of the eigenvalues at F = 0, while the 
open circles indicate those in the limit F —^ oo; the arrows indicate how the eigenvalues evolve as F is increased from 0 to oo. 
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reaches the upper band edge; we can find the precise value of ei where this occurs from the condition 

P(A = -l;£i,rRjc(£i)) =0, (19) 

which yields Si = 1/2. At this precise point, one of the EP2As also touches the band edge and overlaps with the RIC. 
Then for ei > 1/2 the RIC exits the continuum and we find that it splits into a bound state and a virtual bound 
state as shown in Fig. |^e) and (f) for the case £i = 0.6. We also show the evolution of the wave numbers in the 
complex k plane in Fig. [^as the system evolves from £i = 0 to ei = 1.5. Here both move rightward on the real 
axis, excepting that splits into a bound state/virtual bound state pair beyond ei = 1 / 2 . 

Another key difference in the £i 7 ^ 0 case (as well as the eg ^ 0 case) is the appearance of one or (at most) two 
bound states. The bound state properties are discussed in greater detail in Sec. |IV| below. 


IV. FORMAL PROPERTIES OF THE BOUND STATES 


We now turn to a closer investigation of the bound states that appear for various parameter ranges of our PT- 
symmetric prototype model. Our focus here is on the traditional bound states with real energy eigenvalues; however, 
we will briefly comment on the virtual bound states and the localized states with complex eigenvalue at points for 
which they are also relevant to our discussion. 

We first briefly discuss in Sec. |IV A| the parameter ranges for which bound states exist in our prototype model for the 
general case £9 7 ^ 0 , £1 7 ^ 0 and comment on the easiest method for finding bound states for a given set of parameter 
values. Then in Sec. |IVB| we explore the symmetry properties of the wave function for the bound state solutions and 
verify that they satisfy T^T-symmetric boundary conditions; indeed, the virtual bound states (anti-bound states) are 
also T^T-symmetric. Finally, in Sec. IV C we investigate the CPT~ norm [45] for the bound states. 


A. Existence of bound states for the general case £0 7 ^ 0 and £1 7 ^ 0 

As we illustrated in Fig. the bound states of Hermitian tight-binding models can only exist on the real axis of the 
first Riemann sheet, below and above the energy band (this is also true for our present model). For our non-Hermitian 
system, such bound states do not appear in the particular case of cq = ci = 0 as seen in Fig. [^e) and (f); however, 
these do appear for £9 7 ^ 0 or £1 7 ^ 0. For example, one bound state appears above the upper band edge in the case 
£1 > 0 (within a specific range of T values); this state is evidenced in Fig. [^d) by the portion of the trajectory that 
lies on the positive side (Im /c > 0) of the line Re fc = +tt. 

In general, we can write the wave number for the bound states as kj = iKj + Sjir, with Hj > 0 and in which Sj = 0 
for bound states below the lower band edge and 6 j = 1 for bound states above the upper band edge (the same formula 
also holds true for virtual bound states, except that kj < 0 in that case). Then, for a given set of parameter values, 
we can test for the presence of bound states within the spectrum by plugging Xj = = ±e“''^ into P{Xj) = 0 
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Re A Im A 




(a) (b) 



Figure 6: (a) Real part and (b) imaginary part of the roots A of the polynomial -P(A) in Eq.( |ll[ | for ei = —1.1 and eo = 0.05. 
The two real roots 0 < A < 1 become complex at F = 0.45. (c) a region of unbroken PT-symmetry (with all four solutions 
real-valued) in the parameter space (£o,ei,r). 


from Eq, 
solution 


• (CT; ai 
with — 1 


any real solution that yields 0 < Aj < 1 represents a bound state below the lower band edge, and any 
— 1 < A^- < 0 represents a bound state above the upper band edge. Meanwhile, real solutions satisfying 


\j > 1 (Aj < —1) represent virtual bound states below (above) the lower (upper) band edge. In Fig. p3[a, b) we plot 


numerical solutions of (111 for the representative case eg = 0.01 and ei = —1.1. Here we find that there exist two 
bound states below the lower band edge in the parameter domain E < 0.45. There are are also a resonance and an 
anti-resonance in this domain with real part of the energy eigenvalue above the upper band edge. 

We define the unbroken VT-symmetry region as any portion of the parameter space for which all of the solutions 
are real-valued (any combination of bound states and virtual bound states). For example, given £i < 0 we show in 
Fig.§ c) the range of parameter values that yield real values for all four solutions of the dispersion equation. In the 
following Sec. |IVB[ we explicitly demonstrate that both bound states and virtual bound states satisfy PT-symmetric 
boundary conditions. 


B. Verification that real-valued bound states satisfy PT-symmetric boundary conditions 

Here we verify that the real-valued bound states discussed in Sec. |IV A| automatically satisfy PT-symmetric bound¬ 
ary conditions. To accomplish this, we again write the wave number of an arbitrary bound state in the form 
kj = iKj + SjTT, where Sj = 0 for a bound state below the lower band edge and 6 j = 1 for a bound state above 
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the upper band edge. With this formalism, the wave 


equation ([^ for the bound states takes the form 


( for X < - 1 , 

ipj {x) = < V'o for a; = 0 , 

I for a; > 1 . 


( 20 ) 


In order for ipj to be a T^T-symmetric eigenstate of our Hamiltonian H, it must satisfy the condition VTipj = 

Note that at any point we could introduce the state 4’j{x) = e^^^'^ipj{x), which is then an eigenstate of VT with 
eigenvalue 1. Applying the VT operator to the bound-state wave function we obtain 

{ (J*^KX + iS,7,X ^ Q^KX-iS,^X ^ 

■00 for X = 0 , ( 21 ) 

B*e-'^x-iS,nx ^ (J^Kx+^s,7Tx 

where in the last step we have taken advantage of the fact that —tt is physically equivalent to tt in the Brilliuon zone 
structure of our model. If we assume 0 q = FipQ, then we can write the quantity C* jB = B*/C = F as a phase factor 
F = e*®. 

Now, for a T^T^^mmetric solution of our Hamiltonian, we see that we must augment the outgoing boundary 
condition in Eq. ® with an additional condition B = e~'‘^C*, which gives 0(—1) = e“*®'0(l)*. We apply this 
condition to re-write the matrix form of the Schrodinger equation in Eq. as 

/-A + ei+*r 0 \ /0(1)*\ 

-e-*« £0 -1 V'(O) = E{X) 0(0) . (22) 

V 0 -1 -X + ei-iTj VV'(l)/ V0(1)/ 


Taking the determinant of this modified equation yields the exact same condition for discrete eigenvalues P{Xj) = 0 
as we previously encountered at the beginning of Sec. |HI[ Hence, any real-valued bound state of the Hamiltonian F[ 
in Eq. ([^ is automatically an eigenstate of the T^T-symmetry operator with eigenvalue e*®. 

We may obtain the explicit form for the coefficient B = e^^C* from the first and third lines of Eq. (22). For 
simplicity here, let us choose 0 = 0, such that B = C*. We then find the real and imaginary parts of H = Bji -|- iBf 


B„ = 


A (1 -l- Aei) 

1 -k r^A 2 T 2 i^A + e?A 2 


00 ; 


Br = 


TX 


l-kr2A2-k2eiA-ke2A2 


00 , 


(23) 


_ g —/t+2(5j7r 


with A = 

As a final comment, we note that according to the argument we have presented here the wave function for the virtual 
bound states (residing in the second Riemann sheet in the complex energy plane) must also satisfy VT symmetry. 
This can immediately be seen by simply replacing the form of the wave number for the bound state kj = iKj + Sjtt 
by that for the virtual bound states kj = —iKj + Sjir (with kj > 0 in either case) and proceeding with the argument 
as presented above. 

However, we note that the bound states with complex energies are not T^T-symmetric, which can be seen by writing 
the wave number for these states in the form kj = Kj + ilij with H^- < | 7 r| and noting that we can no longer make the 


sign replacement in the final line of Eq. (21 1 as these states reside within the Brillioun zone, rather than at the edges 


as do the bound states and virtual bound states. 


C. CVT norm of the bound states 


We now set out to write the appropriate normalization condition for the bound-state wave function that we obtained 
in Sec. IVB For a non-symmetric Hamiltonian F[, the completeness relation among its eigenstates 0n(a;) assumes 
the form 


= Sn,Tn 


where tpnix) are right eigenstates and ip^{x) are left eigenstates. If H is T^T-symmetric, we identify 


^u{x) =eVT'il^^{x), 


( 24 ) 
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Figure 7: Real part of the roots A of the polynomial Eq.( 11 \ and normalization constant N corresponding to two bound states 
below the lower band edge in the domain 0 < F < 0.45 for ei = —1.1 and eo = 0.05. One bound state has a positive norm and 
the other has a negative norm; the eigenvalues of the C operator are therefore +1 and —1, respectively. A similar picture holds 
for bound states appearing above the upper band edge. 


where the operator C |45] satisfies, in the unbroken region, the three algebraic equations 

[C,VT]=0, [C,H]=0, C^ = l. 

The completeness relation for a T^T-symmetric Hamiltonians then reads 

OO 

[CVT'lpn{x)]llJm{x) = Sn,m- 

X — — 00 


(25) 


(26) 


Since the C operator commutes with the Hamiltonian H, the bound states '4’j{x) of H in the unbroken region must 
also be eigenstates of C and, because = 1, the resulting eigenvalues must be Cj = ±1 |45LI46| . (How the C operator 
might act on the complex-valued solutions in the broken region is a question presently under investigation). 

In order to assign the correct eigenvalue Cj to each eigenstate '4’j{x), we first evaluate the so-called VT norm as 


[VTil)j{x)]il)k{x) = tpj{-x)*i:k{x) = {-lySj^k- 


(27) 


X— — 00 


X — — 00 


We see that the PT-norm is not positive definite ma, with alternating signs ±1 among the bound states '^^(a;); 
hence, we assign the eigenvalues Cn to be ±1 according to the sign of (271 in order to obtain the positive norm 
introduced in Eq. (261. 


In either case, we may write the VT-norm for our bound states given in Eq. (20) as 

OO OO 


N. 


^ — 2K-iX 


A 


(28) 


X — — 00 


x—1 


where Hj is the imaginary component of the wave number from kj = iKj + Sjir. For convenience, we introduce 
B = B/ijjQ', we then obtain 




2^tlR + l 


(29) 


The explicit form of the coefficient B is given by Eq. (23), where A = e~'^ for the bound states below the lower band 
edge and A = —e~^ for the bound states above the upper band edge. In Fig. |^we show the T^T-norm for the two 
bound states previously shown in Figs. |^a,b) that appear below the lower band edge; we see that the VT norm for 
one of these states is positive, giving the eigenvalue of the C operator as 1, while for the other the norm is negative, 
giving the eigenvalue of the C operator as —1. 


V. SCATTERING STATES, THE RESONANCE IN CONTINUUM (RIC), AND PERFECT 

TRANSMISSION 


In this section we consider typical (7^T-asymmetric) scattering boundary conditions for our T^T-symmetric open 
quantum system; we will consider T^T-symmetric scattering solutions later in Sec. |VI[ In Sec. |V A| below we write the 
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generic wave function for the states and then obtain a matrix equation for the relevant scattering coefficients. We 
then explicitly write the transmission and reflection for the simplest case with a pure imaginary potential, namely 
^0 = ^1 = 0- We perform these calculations for both left-to-right and right-to-left scattering and we verify that the 
transmission and reflectance satisfy established relations that generally hold in VT systems. 

In Sec. |VB| we discuss the RIC in detail from the perspective of the scattering wave solutions; here we note that the 
RIC automatically satisfies the Siegert boundary condition for outgoing waves and argue that these states represent 
a resonance between the background continuum and the T^T-symmetric defect potential. 

While the RIC is a discrete state embedded in the scattering continuum, we discuss in Sec. |V C| outstanding 
scattering states, namely perfectly transmissive states, meaning that the transmission is unity while the reflectance 
vanishes. For the simplest case in which the defect potential is pure imaginary with Eq = ei = 0 we demonstrate 
that by appropriately choosing F we can obtain perfect transmission at any given value of k in the spectrum; this 
property approximately holds in the case in which eg and ei are nonzero but take on small values. We also demonstrate 
that the appearance of the perfect transmission state at the band edges coincides with a delocalization transition, 
an observation which may be useful from an engineering perspective as one aims to construct devices with specific 
transmission properties. 

We also examine the case in which not only the transmission is unity, but also there is no phase shift as the scattering 
wave passes through the defect region. In this case the signal is transmitted not only perfectly, but invisibly. 


A. 'PT-asymmetric scattering states 

Here we And scattering states of our T^T-symmetric model Q under the potential ([^. We will limit most of 
the detailed calculations to the simplest case Sg = Si = 0 but the generalization to the case Sg ^ 0 and Si ^ 0 
is straightforward. We first solve the Schrodinger equations for left-to-right scattering by assuming a wave 

function of the form 


Ae^kx 

for X < — 1, 


V’(O) 

for X = 0, 

(30) 

C^ikx 

for X > 1, 



where k resides within the scattering continuum 0 < fc < tt. The term with the coefficient A gives the incoming 
wave, while the B term is the reflected wave and the C term is the transmitted wave. Note that its eigenvalue is real: 
E = —2 cos k. 


We have four parameters A, B, C and ^/'(O) to fix under the three conditions given by the Schrodinger equations (j^- 
Q . Substituting the ansatz (30) into them yields 


/ci -l- *F -|- A —A 0 \ / A \ f ^ ^ 

1 - (1 + eoA + A2) A2 ^( 0 ) = -AH A , (31) 

V 0 -1 l + ie,-tT)xJ \ C J \ 0 J 


Let us limit ourselves from this point to the simplest case eg = ei = 0. Although the overall phase of the wave 
function (301 does not affect physical quantities, it turns out that it is easiest to assume B G M.. We can represent 
the coefficients as 


A = B 


C = B 


V'(O) = B 


i sin k — F^e^*^ cos k 
(F -I- 2 sin fc)r cos k 
i sin k 

(F -I- 2 sin fc)r cos k ’ 
1(1 — jFe*^) sin k 
(F -I- 2 sin k)r cos k ’ 


and thereby obtain the transmission and reflection amplitudes as 

C i sin k 

A i sin k — F^e^*^ cos k ’ 
B (F -I- 2 sin k) F cos k 
A isinfc — F^e^*^ cosfc’ 


(32) 

(33) 

(34) 


(35) 


n 


(36) 
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which are followed by the transmission and reflection probabilities as 


Tl^r := \ti\ = 


Rl^r '■= In I = 


sin^ k 


sin^ fc + (r — 2 sin k)(T + 2 sin k)^^ cos^ k ’ 
(r + 2 sin cos^ k 
sin^ fc + (r — 2 sin fc)(r + 2 sin k)^^ cos^ fc 


(37) 

(38) 


Note that Tr^r + Rr^r is, in general, not unity because we have a source and a sink and therefore the particle 
number is not conserved. Instead, the usual probability conservation relation is replaced by a generalized rule for 
T^T-symmetric systems that relates the left-to-right and right-to-left transmission properties nm, as shown below. 
Hence, we next consider the right-to-left scattering solution given by the ansatz 


ip{x) = 


in which we again have 0 < fc < tt, the D term is the incoming wave, the C term is the reflected wave and the B term 
is the transmission wave. Note again that its eigenvalue is real: E = —2cosfc. 

Again substituting this ansatz into the Schrodinger equations , we obtain 

f B\ / 0 

V'(O) = -AC I A 



for a; < —1, 



for a; = 0, 

(39) 

Cf,ikx D(.-ikx 

for a; > 1, 



H-k(ei-kzr)A -1 0 

A 2 _(i + eoA + A 2 ) 1 

0 —A El — iF -f A, 


\ D ) 


. 1A {si — ^r) A j 


After assuming C G 


this time, we obtain for the case eg = Si = 0 the coefficients as 

isinfc — cosfc 

Jy = (_/ 


B = C 


m = c 


and the amplitudes 






B 

D 

C 

D 


(r- 

- 2 sin fc)F cos fc ’ 


i sin fc 

(r- 

2 sin fc)F cos fc ’ 

*(i 

— zFe®^) sinfc 

(r- 

2 sin fc)F cos fc 


i sin fc 

i sin fc — F^e^®^ cos fc 

. (r 

— 2 sin fc) F cos fc 


i sin fc — cos fc ’ 

which in turn lead to the right-to-left transmission and reflection probabilities 


Tr^r := Id = 


Rr^l ■= Id = 


sin^ fc 


sim fc -I- (r — 2 sin fc)(r + 2 sin fc)r2 cos^ fc 
(r — 2 sin fc)^r^ cos^ fc 
sin^ fc -I- (r — 2 sin fc) (F + 2 sin fc)F2 cos^ fc 


= Tr^r, 


< Rl^r- 


(40) 

(41) 

(42) 

(43) 

(44) 

(45) 

(46) 

(47) 


The left-right asymmetry in (47) comes from the fact that the V and E symmetries are individually broken in our 
system. 

We nonetheless note that we have ti = tj. = t, such that the transmission is equal for the left-to-right and right-to-left 
scattering; this is a general property of VT systems. Further we note that the relations 

t{—k) = t{k) ri{—k) = r^(k) (48) 

are satisfied, which also hold for VT systems in general |1141I115] . Finally, the usual probability conservation property 
for Hermitian systems {T + R = 1) is here replaced by |1141I116] 


|t(fc)| ± |r/(fc)rr(fc)| = 1, 

which we can easily verify by using Eqs. d^, (|M| , ( [44| , and (|4^. This is a result of the fact that Eqs. (|4T|- 


(49) 


are 


obtained from the corresponding Eqs. (32)-(34l by taking the complex conjugate and flipping the sign of fc, which is 
just the VT operation in the wave-number space. Note that the sign choice appearing in Eq. (49) is fixed by the sign 
of the quantity 1 — |t(fc)| = 1 — T; it can easily be shown that for the present case with Eq = £i = 0 the sign changes 
in this quantity generally occur at arcsin(±F/2). 
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(a) Transmission probability 



(b) Reflection probability to the left 



(c) Reflection probability to the right 


Figure 8: (a) The left-to-right transmission probability 1371, which is equal to the right-to-left transmission (461, (b) the left- 
to-left reflection probability (381, and (c) the right-to-right reflection probability (471, all in the simplest case of eq = ei = 0. 
Note that the scale of the vertical axis varies from panel to panel. 


B. Scattering wave perspective on the resonance in continuum (RIC) 


We explain the resonance in continuum (RIC), which we introduced in Sec. Ill as a discrete resonant eigenstate, 
here from the perspective of the scattering states presented in Sec. |V A[ More specifically, we show that the RICs 
appear as singularities in the transmission and reflection probabilities (371, (38), (46), and (47) on the real axis. In 
this sense, it is a discrete state embedded in the scattering continuum. 


Let us recall that in Sec. |III| the RICs appear as the points where two solutions meet on the real energy axis; in 
Fig-i e) for the simplest case Eq = £i = 0i this happens at if = ±-\/2 for F = 1. This is not an exceptional point, 
however, because each of them has a different (real) value of fc; namely, k = ±7r/4, ±37r/4 in the simplest case as is 
shown in Fig. He) or (f)- 

We here show that these points indeed have the properties of the resonant states in the sense that they have 
divergent transmission and reflection probabilities. We plot in Fig. the transmission and reflection probabilities 
both for Eqs. (30) and (39) in the case Sq = Si = 0 for positive 0 < F < 5; they are symmetric with respect to 


k = 0. All probabilities have poles at fc = ±7r/4 and k = ±37r/4 with F = 1, namely for the RICs. These are the only 
instances at which any of the probabilities diverges for real k. 

From this perspective, the poles are indeed discrete states embedded in the scattering continuum. Let us explain 
why they are resonant states in eontinuum. These poles are associated with the zeros of the coefficient A in the wave 
function (30) and the coefficient D in the wave function ([3^, as we can see in Eqs. (35), (36), (44), and ([4^. Exactly 


_ n (f3^ , _ _ _ _ 

at these zeros A = 0 and D — 0, the wave functions (30|) and (|39|) have only outgoing waves, which indeed matches 
the Siegert boundary condition ([^ for resonant states ItellMlTllTgHSS]. Therefore, the poles of the transmission 
and reflectance shown in Fig. are resonance poles. In a resonance state with Refc > 0, particles are ejected from the 
central area and flow away towards x = ±oo; in a state with Re fc < 0, which is historically called an anti-resonance 
state, particles flow into the central area and vanish. We indeed saw this from another perspective in Eq. (17), in 
which the RIC took the form of an outgoing plane wave outside of the central scattering region. 

In the Hermitian scattering problem, the particle number is conserved and hence the Siegert boundary condition 
can only be satisfied at discrete complex values of k and E, which give the location of the resonance poles; it can never 
be satisfied for real values of k and E in the Hermitian case. Under the Siegert boundary condition, the particles flow 
away from the scattering region and hence the particle number in this area decays in time, which can be described only 
by a complex energy eigenvalue 137]. Hence, a resonance pole in the Hermitian case is strongly tied to an eigenstate 
with complex E and k. 

In the UT-symmetric case, however, the particle number is not conserved because we have a source and a sink. 
Hence it is possible for particles to emerge out of the scattering region (or vanish into it) as a stationary state without 
changing the particle number in this region. This is indeed what happens with the resonant states in continuum: 
because Im E = 0 and Im A: = 0 for these poles, the wave function stretches over all space as a stationary state. 

In this sense, it is a remarkable characteristic of UT-symmetric systems that we have resonances resting within 
the real energy continuum. This is why we specifically refer to these as resonance states in continuum; these states 
represent a resonance between the open environment associated with the leads and the UT-symmetric potential. We 
further note that in the special case ei = 0 and F = 1, one or two RICs appear for any values of Eq, although again 
one or the other RIC exits the continuum at Eq = ±1 and splits into a bound state and a virtual bound state, as 
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previously discussed in Sec. 


C. Perfect transmission, invisibility and applications 


In the present subsection, we turn our attention from the RIC poles to the scattering continuum itself. More 
specifically, we now study the system parameters that give rise to perfect transmission such that T = 1 while the 
reflectance vanishes R = 0. We also examine the phase associated with the perfect transmission in order to observe the 
case in which invisibility is obtained, and we comment on several points below that may be useful from an engineering 
perspective. 

The condition to obtain perfect transmission in the left-to-right (right-to-left) case is immediately apparent in 
Eq. ( [30| (Eq. (39l), namely B = 0 {C = 0). This condition is realized whenever the determinant of the 3x3 matrix 
on the left-hand side of Eq. (311 (Eq. @) vanishes. Hence, we obtain left-to-right (right-to-left) perfect transmission 
at a given value k = kxj within the continuum whenever the condition = 0 (Ml(Alj) = 0 ) is satisfied, 

where Axy = and 

-^L,r(A) = (£i T ^r) A^ -f (s^ -|- T^ -f Eq (ci T *r)) A^ -|- £0 (l -|- El -|- A^ -|- (ei -|- T^ -1- Eg (ci i iE)) A -|- £i ± zT = 0. 

(50) 

Since Mx(A) is a quartic polynomial for either case X = L,R, for a given set of parameter values we generally obtain 
four values ki^j for left-to-right and four values k^j for right-to-left perfect transmission. However, some of these 
solutions might turn out to be complex-valued and hence must be discarded. 

In the simplest case Eg = £i = 0, Eq. (501 for the left-to-right transmission gives the factorized form Ml (A) = 
ir(l -I- IA)(1 — iA)(l — ITA — A^). The two linear factors give two T-independent solutions as fcL,i = 7r/2 and 
fcL ,2 = — 7 r /2 with energy E = 0 appearing directly at the center of the transmission band. Meanwhile the quadratic 
factor gives the solutions 


r _ / cos-1 (±^4_r2/2) 

Ld3,4} ^^^-1 (^ 4 -.y /4 _ r2/2) — TT 


cos 


for -2 < T < 0 , 
for 0 < T < 2 . 


(51) 


These four solutions are plotted as the full curves in Fig. |^a); note that the real parts of the wave numbers for the 
eigenvalues Re kj are also plotted as the dotted curves, for reasons that will be described below. From an applications 
perspective, we note for this case we can obtain perfect transmission at any given value of fc S [— 7r,7r] by choosing 
the appropriate value of T. 

So far we have only considered the intensity of the transmitted signal, but not the associated phase information. 
While perfect transmission does maintain the signal intensity, if there is a phase shift between two leads an observer 
may still be able to detect the presence of impurities in the scattering region by performing a time-of-flight mea¬ 
surement m- We can investigate the phase shift for the perfect transmission states by calculating the eigenvector 
(H,'0(0), C')"!' on the left-hand side of the singular matrix (311 for our four perfect transmission solutions. For the 


case of the two perfect transmission values fcL.i = 7r/2 (Al,i = *) and fcL ,2 = —7r/2 (Al ,2 = ~*) we have 



(52) 


(53) 


which yields A = C without any phase shift between the leads and 

0 ( 0 ) = (i±r)H. 

Left-to-right invisibility holds for these two cases; indeed, as we discuss below right-to-left invisibility will also hold 
for k = ± 7 r/ 2 . As a further point, notice from Eq. (531 that we can arbitrarily adjust 0(0) by tuning the parameter 
r and still maintain perfect transmission; for the special case T = =pl we can even eliminate it entirely. 


Performing a similar calculation for the T-dependent perfect transmission states A:l ,3 and /cl ,4 reported in Eq. (511 
reveals that a phase shift is always present for these states, apart from the special case k = ±7r/2 for a specific value 
of T. To summarize, invisibility can only be achieved at fc = ± 7 r /2 for arbitrary T. 


For right-to-left transmission in the simplest case Eg = £i = 0, we obtain from the lower sign in Eq. (50 1 the 
factorized equation Mr(A) = —ir(l -|- iA)(l — zA)(l -I- iTA — A^). Hence for the right-to-left transmission we again 
obtain perfect transmission at = 7 r /2 and = — 7 >'/ 2 , but for the T-dependent perfect transmission states we 
have instead 


«R,{3,4} 


cos 


-1 


(±V4-r2/2) - TT 


cos-1 (±^4_r2/2) 


for -2 < T < 0 , 
for 0 < T < 2 , 


(54) 
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(c) 


(d) 


Figure 9: The wave numbers for perfect transmission (full curves) as a function of F for fixed values of so and ei; the wave 
numbers Re kj for discrete eigenvalues are also shown in the background with dotted curves, except for bound states which 
appear as full lines at either k = 0 (lower band edge) or ±7r (upper band edge). We used the following parameter values: 
(a) left-to-right transmission for eo = £i = 0, (b) right-to-left transmission for eo = £i = 0, (c) left-to-right transmission for 
So = 0, £i = 0.08, and (d) left-to-right transmission for so = —0.1, £i = 0. 


as plotted in Fig. ib)- Notice that the expressions ( [54| are just reversed from the left-to-right case reported in 
Eq. (511; this is a natural result for a T^T-symmetric system since switching our scattering orientation amounts 
to swapping the position of the gain and loss impurities. Again in this case we can evaluate the scattering wave 
coefficients to find that D = B and hence we have no phase shift for the perfect transmission states fcR,^{i, 2 } = i7r/2; 
this shows that these states yield bi-directional invisibility. 

Here the response at the site 0 is given by 


^(0) = (lTr)f9. (55) 

If we choose, say F = —1, the invisible right-to-left wave has a finite value of '0(0) at fc = -l-7r/2, but from Eq. ( [5^ 
the left-to-right transmission gives no response at the same frequency. Hence the site 0 in this scenario can act as a 
kind of switch that responds to the otherwise invisible signal transmission in one direction but ignores signals from 
the other direction. 

In Fig.j^c) and (d) we plot the left-to-right perfect transmission for more general parameter values with Eq = 0, Ei = 
0.08 and Eq = —0.1, Ei = 0, respectively. We note that the range of the continuum that is capable of supporting perfect 
transmission has been reduced slightly in comparison to the ‘cleaner’ case in Figs. |^a) and (b) for these relatively 
small values of the impurity energies. For larger impurity values the range of coverage for perfect transmission is 
further reduced. 

We also note the following connection between the perfect transmission scattering states and the bound states of 
the discrete spectrum. Notice that in the background of Fig. [^a)-(d) we have plotted the real part of the wave 
numbers Refc^- for the discrete eigenvalues as the dotted curves. However, the wave numbers for the bound states are 
marked with full lines appearing at either fc = 0 (lower continuum edge) or k = ±7r (upper continuum edge). Focusing 
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on Fig. I^c) and (d), we note that the appearance of the perfect transmission scattering state exactly coincides with 
the appearance (or disappearance) of a bound state at the edges of the Brillouin zone k = 0 or k = ±7r; see the 
footnote |120| . 

This is a rather intuitive result if we consider the behavior of the bound-state wave function at the delocalization 
transition, where it brushes against one of the band edges before becoming a virtual bound state in the second 
Riemann sheet. On one side of this transition we have a bound state with a wave function that is localized in the 
defect region, while on the other side we have a virtual bound state with a wave function that diverges into the leads 
(one can even view this state as being localized in the leads [33]). At precisely the transition between these two states, 
we have a wave function that spreads out evenly throughout the chain, which here supports perfect transmission from 
one lead to the other. 

This explains the connection between the delocalization transition and the appearance of the perfect transmission 
state at either edge of the scattering continuum, which may provide an intuitive approach to designing systems with 
desired transport properties. We notice a somewhat similar transition appears in Ref. |108| . 


VI. RT-SYMMETRIC SCATTERING WAVE SOLUTIONS 

This section is devoted to more mathematical interest. We here investigate the T^T-symmetric properties of the 
scattering wave solutions. Previously we demonstrated in Sec. |III| and Sec. |IV| that the discrete states satisfy VT- 
symmetric boundary conditions in certain regions of parameter space; however, we note that despite the fact that the 
eigenvalue E = —2cosfc associated with the scattering states is always real, in Sec. these states generically satisfied 
PT-asymmetric boundary conditions. 

This motivates us to investigate whether or not the scattering states themselves can obey PT-symmetric boundary 
conditions. We will show in Secs. |VI A| and |VI B| that we indeed always have PT-symmetric scattering states. In 
the former, we will PT-symmetrize the scattering wave function previously obtained in Sec. ^ In the latter, we 
will present a more direct and systematic way of finding a PT-symmetric scattering wave with the use of the Jost 
solutions. We will then introduce the concept of the VT current in Sec. |VI C[ 


A. 'PT-symmetrization of the scattering wave solutions 


We can construct a PT-symmetric solution out of an asymmetric solution ijj{x) by writing 

ip-prix) = ijj{x) + VTi^ix). 


(56) 


Let us apply this strategy to the scattering solution (30l. The VT transformation results in the changes i ^ —i and 
X —>■ —X as well as the complex conjugation of coefficients A and C; we recall that B is real by assumption in the 
relations (32|-(34) for the simplest case Eq = ei = 0. Note that the real-valued wave number k (the momentum) is 


invariant under the action of VT, since both V and T result in a sign flip separately for this quantity. Substituting 
our result into Eq. (561 we obtain a VT -symmetric solution for 0 < k < tt given by 


V'pr 


(A -h ■ 

(x) = { + jP{0 )* 


Be 


— ika. 


with the relations (32 h( 34| producing 


{A* + + Be 


A + C* = -B 


r, — ika. 


for X < — 1, 
for X = 0, 
for X > 1 


o2zfcT 


V’(o) -f V'(O)* =B 


2 sin k + T’ 
2 sin k 


2 sin fc -f r 

(L), 


(57) 

(58) 

(59) 


for the simplest case Eq = ei = 0. Note that the component '>p}p^{0) is real. If we choose the normalization as 
= 1 we obtain 




Ak(x+2) 


2 sin k 


+ 1 + 


Ak(x-2) 


2 sinfc 


+ 1 + 


2 sin k 


2 sinfc 


^ — ikx 


^—ikx 


for X < —1, 
for X = 0, 
for X > 1 


( 60 ) 
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for 0 < fc < TT, as our first T^T-symmetric solution. 

We can instead start from the right-to-left scattering wave Eq. (39l. Following a similar procedure as above, we 
obtain 


+ {B + 

= “I V'(O) + ^(0)* 


for X < —1, 
for a; = 0 , 
for X > 1 


for 0 < fc < TT, where we used the assumption C G K in this case. In the simplest case £o = £i = 0 we have 

-2ifcp 

B + D* = C- 


^{0)+i;{0)* = C 


(¥i) 

After again choosing our normalization such that (f>^^{0) = 1 we obtain 


2 sin A: — r 
2 sin k 
2 sin k — T 


(61) 

(62) 

(63) 


= < 


1 - 


1 - 


2 sin k 

r 

2 sin k 


Akx 


ikx 


^-ik{x+2) 


2 sin k 


^-ik(x-2) 


2 sin k 


for X < — 1 , 
for X = 0 , 
for X > 1 


(64) 


for 0 < A: < TT, as our second T^T-symmetric solution, which is indeed obtained simply by flipping the sign of k in the 
first solution (60l: We therefore conclude that the solution (601 holds for the entire first Brillouin zone —tt < A: < tt. 


B. Jost solutions 


The solutions in the previous subsection |VI A| seem somewhat strange because of the asymmetry with respect to 
the inversion of k. In this subsection we obtain an alternative T^T-symmetric solution by directly finding the Jost 
solutions of the original Schrodinger equation We will find a solution of a more symmetric form, which is 

indeed a superposition of the solutions (60) and ( 64]). We again restrict ourselves to the simplest case Eq = Si = 0. 

Let us briefly overview the way of constructing a scattering wave solution out of the Jost solutions. When the 
potential vanishes far away from the origin, we can assume a solution of the form of a plane wave there. The solutions 
thus defined under the boundary conditions |118l I119| 


/±(x) = ae 


^ikx 


as X 


oo 


(65) 


with an appropriate constant a are called the Jost solutions. They, however, do not generally satisfy boundary 
conditions at the origin. We therefore take a superposition of the two Jost solutions so that it may satisfy the 
boundary conditions at the origin, which yields a scattering wave solution. 


Since the potential vanishes for x > 2 in the present case, we can use Eq. (65) in the region x > 1. Let us now 
construct a T^T-symmetric Jost solutions. Since we set 


/±(x) = a* e 


* ^^ikx 


as x —)• —cx), 


( 66 ) 


which we can use in the region x < — 1. These Jost solutions, however, do not satisfy the boundary conditions at 
X = 0. Indeed, the Schrodinger equation for x = 1, namely Eq. (0, gives 


/±( 0 ) = -/±( 2 ) + (-E(fc)-*r)/±(l) 


= a|-e=^2i/c 


= a(l - zTe^*'') 

for El = 0, while the Schrodinger equation for x = —1, Eq. ([^, gives 

/±(0) = -/±(-2) + {-E{k) + *r)/±(-l) 

= a*[—e 


(67) 




( 68 ) 
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We can make Eqs. (671 and (681 continuous at the origin by choosing 


a= l + 


(69) 


which makes /±(0) = 1 ± 2r sinfc + F^, but the resulting solution f±{x) does not satisfy the Schrodinger equation for 
a; = 0, Eq. (|^ (even after setting £o = 0 for the present case). 

The physical solution that satisfies the Schrodinger equation ([^ must be given by a linear combination of the two 
Jost solutions: 


(pvrix) = A+f+{x) + A_f_{x) (70) 

with two superposing coefficients A±, which we set so that (j)'p'f{x) may satisfy Eq. (§. Let us define the Jost functions 
(not to be confused with the Jost solutions) by 

F±{k) = -/±(-l) - /±(1) - E(fc)/±(0) 

= -(1 - - (1 + + (e*'^ + e-*'=)(l - fre=^*'=)(l + 

= 2r(r±2sin/c)cosfc. (71) 

The Schrodinger equation ([^ is then reduced to 

A+F+{k) + A_F_{k) = 0, (72) 


which fixes the ratio between A+ and A-. 

By normalizing the function by (j)-p'f{0) = 1, we can express the final result as follows: 


(fvrix) = < 


1 - 


1 - iTe 


ik 


1 - 


2 sin k ] 2 

r \ l+iTe"*'' 


+ ( 1 


r 


2 sin A: 




2 sin k 

r 


—ik 


l-iTe 
-( 

2 

1 + 


2 sin k 


^ — ikx 


—ikx 


This is indeed a linear combination of Eqs. (60) and (|64|) in the form 

r 




) A'‘F) + 2 (' 2sint) FFy. 


2 sin k, 

however, the domain extends over the entire first Brillouin zone —tt < fc < tt. 


for x < —1, 
for a: = 0, 
for a: > 1. 


(73) 


(74) 


C. ■PT-current 

Because we have a source and a sink, the discrete states generally do not conserve the particle number and the 
scattering states do not conserve the current. We here, however, introduce a current that is conserved for a VF- 
symmetric scattering state, which we refer to as the T^T-current. 

The standard current is defined in a one-dimensional continuous space as 

j = Re (il}[xYp%f{x)) = j. (^fj{x)* ^-ipix) - 'ipix)^'ip{x)*^ , (75) 

which would normally be independent of x, but this does not generally hold true in a T^T-symmetric non-Hermitian 
system. We here instead introduce the T^T-current 

jvT = ^ 

We can prove that the T^T-current is independent of x for an eigenfunction fifx) with real eigenvalue E of the 
Hamiltonian 


H-pT = — 




( 77 ) 
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with VTVp'f{x) = V-pri—x)* = Vp-p{x). Computing the x derivative of the T^T-current (76), we indeed have 


dx'' 


2 i 


1Z^vt{x) = - ( i^{x)* ^i^{-x) - ^{-x)^i^{x) 


dx"^ 


= ^ [^{xY (Vpri-x) - E) Y[-x) - Y{-x) (VvrixY - E) Y{xY] = 0. 


(78) 


Notice that it vanishes identically for a T^T-symmetric eigenfunction, because we then have YixY — Y{~x) in Eq. 
In the discretized space of the one-dimensional tight-binding model, the standard current is given by 


j = ^ [Y{xY {Y{x +1) - Yix)) - Y{x) {Y{x + 1)* - Y{xY)] 

= (YixYYix +1) - Y{x)Y{x + 1)*), (79) 


while the T^T-current is given by 


jvT = ^ {Y{xYY{-x - 1 ) - Y{-x)Y{x + 1 )*) • 


For a T^T-asymmetric left-to-right scattering state of the form (30), the (traditional) current (79) is 

j = sin k X 


\A\^ + \B\^ 


for X < — 2 , 
for cc > 1 , 


(80) 


(81) 


which are generally not equal along the two leads as we showed in Sec.|^ The T^T-current (80), on the other hand, is 


jpj- = sin k X 



(r -I- 2sinA;)r cosfc’ 


for X < — 2 , 
for X > 1 , 


(82) 

(83) 


which is conserved on both sides of the scattering region. We plot the T^T-current in Fig. [T0| the singularities 
appearing here correspond to the left-to-right perfect transmission states previously shown in Fig. |!^a). 
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VII. CONCLUSION 


In this paper we have combined two types of non-Hermitian systems, open quantum systems and T^T-symmetric 
systems, in order to study a simple example of a T^T-symmetric open quantum system. This system took the form of 
a tight-binding model with a T^T-symmetric defect potential as shown in Fig. which might be physically realized as 
an optical lattice array or approximated in a variety of VT systems with a defect scattering center |CT1 [501 1102] . We 
used this model to illustrate a number of quite general features of T^T-symmetric open quantum systems, including 
properties of the discrete spectrum as well as the scattering states. 

In Sec. Ill we studied the resonance state in continuum (RIC) as a feature of the discrete spectrum, illustrating 
that it represents a resonance state appearing directly within the conduction band (scattering continuum) as it crosses 
from the second Riemann sheet of the complex energy plane into the first sheet. In Sec. jlllBj we showed this state 
takes the form of an outgoing plane wave from the impurity region into the leads. As a result, this feature also appears 
in the scattering continuum and hence we returned to it in Sec. fV] in which we studied the scattering properties of 
the system; in Sec. jV Bj we described that the RIC represents a resonance between the open environment associated 
with the leads of the system and the T^T-symmetric defect potential. In this sense, the RIC can be viewed as a quite 
particular feature of T^T-symmetric open quantum systems. 

We also showed in Sec. jlllDj that the RIC may exit the conduction band as we modify the system parameters by 
splitting into a bound state and a virtual bound state at the edge of the continuum; we believe that this effect should 
be experimentally observable, perhaps in a modified version of the experiments presented in Ref. [5^ or Ref. |60| . 
We also point out that it has been previously illustrated that a bound state (or virtual bound state) appearing near 
the edge of the continuum should generally result in an enhancement of the long-time non-exponential decay that is 
known to appear in open quantum systems Since both a bound state and a virtual bound state appear near 

the band edge in the case when the RIC exits the continuum, this may result in an even greater enhancement of the 
non-exponential effect that could also offer an interesting basis for experimental study. 

We illustrated another key difference from Hermitian open quantum systems in that complex-valued solutions are 
allowed to appear in the first sheet of the complex energy plane in the T^T-symmetric case. These appear as pairs 
of localized states, one with an amplifying characteristic and the other with an absorbing characteristic, as observed 
experimentally in Ref. j5n|. We also pointed out that some of these states may behave as quasi-bound for large values 
of the VT defect parameter T; these again might be observable in a system that imitates our defect potential. 

We evaluated general scattering properties of T^T-symmetric open quantum systems in S ec.[V| in which we calculated 
the transmission and reflectance for a T^T-asymmetric scattering wave solution in Sec. jV A] and verified that these 


satisfy known symmetry relations |114H116] . We also studied perfect transmission states in Sec. VC with invisible 
solutions as a subset of these, and illustrated a connection between perfect transmission at the band edges and a 
delocalization transition of the bound state in the discrete spectrum. 

After noting that the eigenvalue E = —2cosfc associated with the scattering states is always real, in Sec. jVIj we used 
our model as a mathematical prototype to illustrate the construction of scattering wave solutions that themselves 
satisfy T^T-symmetric boundary conditions (just as the bound state in the discrete spectrum is well known to satisfy 
such boundary conditions, as we have illustrated in Sec. IVI. In Sec. VIC we wrote the T^T-current for these solutions 


and pointed out that the previously studied perfect transmission states appeared as a special case. 
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Appendix A: EP eigenvalue expansion in the case eq = £i = 0 


Here we briefly detail the eigenvalue expansions obtained in the vicinity of EP2As (Eqs. ((T4|) and (13l) and EP2Bs 
(Eqs. (151 and ( [T^ ) for the case Eq = ei = 0, following a variation on the method developed in Ref. 77]. First we 
find it useful to rewrite the polynomial equation P(A) = 0 from Eq. (11 \ directly in terms of the energy eigenvalue E. 
We accomplish through the substitution X = —{E + \/E‘^ — 4), which yields the equivalent equation p{E:j) = 0, with 


p{E) = r^E* -f (F'^ - 4r^ -1)E^ + 4. 


(Al) 
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(Note that we have chosen different labeling j for the solutions of this alternative form of the dispersion equation in 
order to emphasize that there is no consistent labeling that will hold between the sets of solutions as we cross the 
EP dH]). 

The basic idea for our calculation is that we will take advantage of the fact that the derivative of the eigenvalues 
blow up at the EPs to study the system properties nearby. First we take a full derivative of the polynomial equation 
dp/dP = 0 and re-arrange to obtain 


2 rE 4 -t 4 r (r^ - 2) 

dE/dT 


-f 2E [2r^E^ -k r'‘ - dP^ - 1] = 0. 


(A2) 


Since dE/dV diverges, we obtain the a useful relationship between E = E and P = P at the EP by setting the second 
term on the RHS above to zero, which yields 


E{T = P) = ± 


Vi -h 4r2 - r4 

Tif ■ 


(A3) 


We can then plug this formula back into the original polynomial dispersion given in Eq. (A1 \ to find the locations of 
the EPs in parameter space as P = PPa and P = APb, where 


f A = V 2 - 1 


f B = 1 + V 2 . 


(A4) 


We then find the locations of the eigenvalue coalescence points by plugging these values back into Eq. (A31 to find 
E{fA) = ± \Ea\, P^(—Pa) = ± I-^aI, ^'(fB) = I^bI, and if(—P b) = ±i |Ab|, with 


^^A = 


^ 2(1 + 72 ) 


Eb = W2 V2-1 . 


Following Ref. we can now write a generic expansion in the vicinity of each the EPs as 


EA(r) = Ea +ai,2Vr2-fi, 


and 


Ea(t) = -|Ea|+ a3.4Vr2-fi; 


£;B(r) = ^IEbI +/3i,2-\/r2-f|, 


EB(r) = Eb +/33.4Vr2-f|. 


(A5) 


(A6) 


(A7) 


To find the expansion coefficients ai^ 2 ) for example, we define = P^ — P^, plug this into Eq. (A2), and expand in 
powers of A. Carrying this out for both cases we obtain 


cri,2 = 0^3 .4 = 


21/V-l + V 2 ’ 


A ,2 = PsA = ±i 


21/4^1 + 72’ 


(AS) 


Puttin g Eqs. ( A5 l and (^1 into Eq. (Ml and Eq. ( |A7[ ), we obtain the expansion associated with the EPs as reported 
in Sec. lmil 


Appendix B: Properties of complex localized states in Region IV for the case £o = £i = 0 


In this appendix we detail the properties of the complex localized states in Region IV (see Fig. |3|)b)) for the case 
£o = = 0; as discussed near the end of Sec. [Ill A| and Sec. [Ill C| of the main text. We can generally expect the 

condition A ^ 1 to hold throughout this region of the parameter space. Hence we begin by expanding the solutions 
of Aj, which are reported in Eq. (121, in powers of l/P to obtain 


A 


1,4 




A 


2.3 


±*|l-f2 


(Bl) 
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We use Ej = —{\j + to obtain expansions for the energy eigenvalues immediately as 


£^1,4 


±i 




E2,3. 



(B2) 


as reported in the main text. 

To elucidate the asymptotic properties of the wave function for these eigenvalues, we make use of the first and third 
rows from Eq. to write 


^(t1) ^ 1 ^ A 

?/'(0) -A±zr-E(A) liiAT' 


(B3) 


First, let us evaluate the localization properties for 'tpi associated with the eigenvalue Ei ~ *r, which appears to be 
the uncoupled gain site in the limit T —^ oo. The calculation for the wave function proceeds along similar lines. 
Applying AiT = i{l + l/T^) we find 


V^i(-l) _ _ -p _ . J_ 

V'i(o) * ’ V'i(o) '^*2r- 

Choosing our normalization such that '0i(—1) = 1, we have '0i(O) « i/T and 
Ai = « i/T to write the wave function ([^ as 


r for a; < -1 

~ N p for a; = 0 , 

I ^ for a; > 1 


or 


(B4) 

— l/(2r^). We then use 


(B5) 




(B6) 


which shows that the state is localized around the gain site x = — 1 with the localization length l/logT. A similar 
calculation for E 4 « —iT shows that the wave function for this eigenvalue is localized around the loss site x = 1. They 
become sharper and sharper as we increase T. 

For the eigenvalues ^2^3 in Region IV, we can use A2,3 « ±J to find 


V>m(Tl) _ -A 

V’m(o) ^ r 


(B7) 


for both eigenvalues with m = 2,3. This indicates that the site 0 is the localization center for the wave functions 
Let us therefore normalize the wave functions according to '0^(0) = 1- Because the wave numbers are expanded as 


we obtain the wave functions as 


fc2,3 =-jlogA 2,3 « -ilog ^ 

TT . 1 


( zp ie±*’^l“l/2e-l“l/r' for x < -1 
02,3(2;) « 1 for X = 0 

^ lg±i7rx/2g-x/r" for X > 1 


(B8) 


(B9) 


which shows that the localization length is r^/2. In other words, they become broader and broader as we increase T. 
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Note that we can explicitly show that the appearance of the perfect transmission state and the bound state delocalization 
do indeed occur at the same point in parameter space in the following manner. On one hand, we can exactly locate the 
delocalization transition at the band edges A = ±1 by plugging into the dispersion polynomial Eq. (Ill as P{\ = ±1) = 0; 
on the other hand we can find the appearance of the perfect transmission states at the band edges by plugging into Eq. 
( |50| ) as Ml(±1) = Mr(± 1) = 0, which can be analytically solved in either case. Doing so yields the exact same points in 
parameter space. 



